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We examine theoretically the role of spin-waves on current-induced domain wall dynamics in a 
ferromagnetic wire. At room temperature, we find that an interaction between the domain wall 
and the spin waves appears when there is a finite difference between the domain wall velocity 
±o and the spin current u. Three important consequences of this interaction are found. Firstly, 
spin-wave emission leads to a Landau-type damping of the current-induced domain wall motion 
towards restoring the solution xq = u, where spin angular momentum is perfectly transfered from 
the conduction electrons to the domain wall. Secondly, the interaction leads to a modification of the 
domain wall width and mass, proportional to the kinetic energy of the domain wall. Thirdly, the 
coupling by the electrical current between the domain wall and the spin waves leads to temperature- 
dependent effective wall mass. 
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I. INTRODUCTION 

The advent of giant magnetoresistance, magnetic tun- 
nel junctions, and spin-transfer torque in magnetic het- 
erostructures has led to proposals of novel applications in 
which magnetic domain walls are manipulated by electri- 
cal currents^ instead of magnetic fields!^ Experimental 
studies based on such concepts have been made possi- 
ble by vast improvements in nanofabrication techniques, 
which allow for more prec ise control over domain wall 
nucleation and propagation While current-driven wall 
motion is an important means of realizing potential ap- 
plications, the threshold current for such motion still re- 
mains prohibitively high for use in integrated circuits P 
As a consequence, strategies are being actively sought to 
simultaneously achieve low threshold current densities in 
combination with high-speed domain wall motion. 

From the point of view of fundamental physics, 
current-driven domain wall motion has attracted much 
interest because it associates a complex spin-dependent 
transport problem with nonlinear magnetization dynam- 
ics. This is equally true for ferromagnets based on 3d 
transition metals, such as iron, nickel, cobalt, and asso- 
ciated alloys, as for dilute magnetic semiconductors such 
as (Ga,Mn)As. From a theoretical perspective, the prob- 
lem lies in computing the correct torques exerted on the 
magnetization by the conduction electron spins. If one 
assumes that the conduction electron spins, propagating 
with an effective drift velocity it, track perfectly the local 
magnetization along their passage through the domain 
wall, one finds an additional torque on the magnetiza- 
tion M of the form 

T a = -(u-V)M, (1) 

which is often referred to as the "adiabatic" contri- 



bution of spin-transfer. This term is well-understood 
and has been reproduced from different transport the- 
or j es [UUMM!]] The magnitude of the effective drift ve- 
locity is given by u — jPg/is /(2eM s ), where j is the 
charge current density, P is the spin polarization, /is is 
the Bohr magneton, e is the electronic charge, and M s 
is the saturation magnetization. An outstanding prob- 
lem of importance concerns the origin of the so-called 
"non-adiabatic" contribution^ 

T na = A M x [(« • V)M] , (2) 

which has been found to be necessary to describe some 
experimental data. M s is the saturation magnetization 
and the dimensionless coefficient /3 characterizes the mag- 
nitude of the non-adiabatic contribution. 

It is possible to gain good insight into the physics 
of current-driven motion through the one-dimensional 
model (ODM) for domain- wa ll dyn amics. This model 
was much studied in the 1970s Poland later adapted to 
the case of exchange to rques du e to coupling to conduc- 
tion electrons by Berger^* 16 * 17 ! and Tatara and Kohno.^ 
The ODM is derived from the Landau-Lifshitz-Gilbert 
(LLG) equation of motion for magnetization dynamics 
with the current-driven terms 

dM „ a „, dM m m 

— = lQ H cS x M + —M x — + T a + T na , 3 
at M s at 

where 70 is the gyromagnetic constant and a is the 
Gilbert damping constant. The relevance of Gilbert 
damping regarding domain wall motion was dicussed 
by Stiles et aZP^ By assuming that the domain wall 
shape remains rigid during propagation, it is possible to 
parametrize the dynamics in terms of only two conju- 
gated coordinates: the domain wall position, xq, and its 
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conjugate momentum, p. By assuming that the external 
forces acting on the wall are sufficiently weak such that 
the wall shape remains rigid, the equations of motion in 
this limit can be written as,^ 

dp 2aK±mdx 2f3K ± m 

dt = — ^^r + ^^ M+/zc + /pin; (4) 

dxo _ aS dp p 

dt ~~ 2K ± mdt m' 1 ' 

where m is the domain wall mass, S is the spin angular 
momentum at each individual magnetic site and K± is 
the transverse anisotropy energy. The external magnetic 
field and a pinning potential, due to intrinsic defects or 
artificial pinning centers, generate the additional forces 
/ze+/pin on the domain wall, respectively. The first equa- 
tion Q relates the domain wall acceleration dp/dt to the 
total force. The second equation (|5| relates the domain 
wall velocity to the domain wall momentum. Neglecting 
damping, one sees that p is related to the relative veloc- 
ity ±o — u through the Doring mass of the domain wall 
m = pj (xq — u). 

The importance of the non-adiabatic "/3-term" is made 
explicit in Eq. [4] where its contribution as an effective 
magnetic field can be immediately seen. It has been 
shown in previous studies that the existence of (3 leads 
to different qualitative dynamics for the wall motionP If 
2/3 > \J H p /H±, with H p being the extrinsic pinning field 
and H±_ the transverse anisotropy field, the domain wall 
depins for u > \"foH p /2[3. Therefore in the weak pinning 
limit, the larger is the (3 term, the smaller is the critical 
current density™ 

However, as we have indicated above, the physical 
origin of this non-adiabatic term is still an open issue 
subject to spirited debate.^ In one line of inquiry, dif- 
ferent authors have sought to associate (3 with the vis- 
cous damping coefficient a, since both parameters de- 
scribe dissipative processes.^ Barnes and Maekawa con- 
tend that (3 and a are equal because of Galilean invari- 
ance|22l while Kohno et aZ.^Duine et aZ.j^and Piechon 
and Thiaville^ have found that [3 and a are not equal 
in general. In a different picture, Tatara and Kohno as- 
sociate (3 with ballistic domain wall resistance, 7 which 
is independent of a and depends only on the transport 
properties of the system. Much of the difficulty in reach- 
ing a consensus is therefore related to the complexity in 
defining the /3-term theoretically and in measuring it ex- 
perimentally 

The present study is motivated by the hypothesis that 
the interaction between the domain-wall and spin-waves 
produces a term similar to the non-adiabatic term, but 
in the presence of only the adiabatic component of spin- 
transfer. The role of spin waves on field-driven domain 
wall dynamics h as been examined by a number of au- 
thors in the pastJ2§12ZI2§l but their role on current-driven 
wall dynamics has not been studied in much detail the- 
oretically. While most theories on the /3-term have fo- 
cused on the transport properties of the conduction elec- 
trons, few studies have considered the motion of the non- 



equilibrium magnetization by taking into account the 
fluctuations. Nevertheless, the interplay between spin 
waves and the domain wall should be important for at 
least two reasons. Firstly, thermal spin waves account 
for a decrease in the magnetization which can be im- 
portant if the system temperature approaches the Curie 
temperature. This is certainly the case in dilute mag- 
netic semiconductors. Secondly, the spin waves act as 
a thermal bath with which energy can be exchanged 
with the domain wall. Indeed, the importance of spin- 
waves as a channel for energy dissipation in magnetic 
system has long been recognized. In the context of the 
ferromagnetic resonance, for example, two-, three- and 
four-magnon processes have been shown to be crucial for 
explaining resonance linewidths of ferromagnetic insula- 
tors!^ In the context of domain wall motion, Bouzidi 
and Suhl have showed that power is diverted from the 
domain wall motion through the amplification of some 
thermal spin wavesP 

Recent experimental studies suggest that current- 
induced domain wall motion may depend strongly on 
the temperature. | 30 | 31 | 32 | 33 | 34 | 35 l Experiments on current- 
induced domain wall motion in metallic devices are gen- 
erally performed at room temperature, but recently sev- 
eral measurements have been reported over a range of 
temperature from several dozens to a few hundreds of 
degrees KelvinP^l Studies on the temperature are likely 
to bring detailed information on the current induced do- 
main wall dynamics. The actual temperature of a fer- 
romagnetic wire along which a charge current flows is 
generally modified by Joule heating, and may vary much 
from one sample to another depending on the efficiency 
by which heat is drained out. As the current density 
required for pushing a domain wall in a ferromagnetic 
metal is usually quite high, j sa 10 12 A/m 2 , the increase 
in the temperature due to Joule heating may even ap- 
proach the Curie temperature Tj^l This could induce 
a drastic effect on the saturated magnetization. Similar 
heating effects may also appe ar in nanowires involving 
magnetic semiconductors. 31 32 Laufcnberger et al. found 
current driven domain wall motion to be less efficient 
by 50 % when temperature increases by 200 KP These 
authors suggest that this loss of efficiency is due to the 
excitation of spin waves. 

In this paper, we study the role of spin waves on 
current-driven domain wall motion by extending the ap- 
proach used by Bouzidi and SuhlJ^El which associates 
some basic i deas for m the theory of solitons^ with spin- 
wave theory. ^ ! 37 ! 38 ! The coupling between the domain 
wall and the thermal bath of the spin waves, w hich o rig- 
inates in the kinetic part of the spin lagrangianj 39 ! 40 ! has 
a number of consequences on the current-induced do- 
main wall motion. It leads to a new dissipation chan- 
nel, whereby magnons can be absorbed or emitted as the 
domain- wall propagates. This dissipation channel relaxes 
the domain wall dynamics towards the solution xq = u, 
where the domain wall velocity and the conduction elec- 
tron spin current are identical and Galilean invariance is 
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FIG. 1: Geometry. The wire is along x axis and the static 
domain wall profile is in the plane (x, y) because of a strong 
perpendicular anisotropy along z. The spherical polar coor- 
dinates are defined with respect to the z direction. 



restored. This dissipation process is somewhat analogous 
to Landau-damping in plasmas^. The coupling between 
the spin waves and the current driven domain wall also 
results in stochastic forces in addition to damping. These 
stochastic forces are weakly correlated at the time scale 
of the domain wall thus behave as a white noise. Besides 
if galilean invariance is not satisfied u ^ x'o, the flow of 
the spin current across the wire will reduce the domain 
wall width and will renormalize the energy of the sys- 
tem. This renormalization of the system energy can be 
re-interpreted as a modification of the domain wall mass, 
which becomes temperature-dependent through the in- 
teraction with the spin- waves. 

This paper is organized as follows. The spin- wave 
cigcnmodes of the domain wall are determined in Section 
II. In Section III the ODM of current induced domain 
wall dynamics is generalized by taking the spin waves 
into account. In Section IV, the damping of domain wall 
motion through radiation of magnons is presented. This 
radiation leads to both a-like and /3-like terms, which 
are both proportional to the domain wall kinetic energy 
p 2 /2m. The change in domain wall width by the electri- 
cal current is calculated in Section V and subsequently 
interpreted as a renormalization of the domain wall mass. 
The response of the spin waves to the domain wall dis- 
placement and to the spin-transfer torque is investigated 
in Section VI. The renormalization of the domain wall 
mass, as a result of this response, is then estimated nu- 
merically. In Section VII, we present some discussion 
and concluding remarks, as well as offering suggestions 
for new experiments that are designed to test the main re- 
sults of our theory. The Green functions used for our cal- 
culations and the integral equation used for determining 
the spin-wave response are presented in the Appendix. 



II. EIGENMODES OF A BLOCH DOMAIN 
WALL 



We consider a ferromagnetic wire lying along the x 
axis, with an axis of easy anisotropy K u along the x 
direction and an axis of hard perpendicular (transverse) 
anisotropy K± along the z direction (see Fig. [I]). The 
orientation of the localized spins is described in spherical 



coordinates within a continuum approximation by means 
of a field m(x,t) — (sin 9 cos <p, sin 9 sin <fi, cos 9), where 
0(x,t) and 4>(x,t) have a space and time dependence. In 
the absence of a conduction electron charge current, the 
magnetic energy H of the system is 

H = J d 3 r {A ((V6>) 2 + sin 2 9(X7cf>) 2 ) 

-K u sin 2 9 cos 2 <j) + K±_ cos 2 9} , (6) 

where A denotes the exchange coupling. In the follow- 
ing we assume K± ^> K u . This condition applies well to 
thin permalloy nanowires, whose in-plane (magnetocrys- 
tallinc) anisotropics are very small compared to the per- 
pendicular (demagnetizing) energy. The anisotropy con- 
stant K u describes the shape anisotropy in the plane of 
the wire and is as small as a few Oersteds, whereas the 
anisotropy constant K± describes the demagnetizing field 
and is about 4irM s « 10 kOe. 

It is stated by the principle of stationary action that 
the energy of the magnetic system at equilibrium is min- 
imized. In other words, 



SH, 



5H r . 



which leads to 



70 = 



.4 



dx 2 



0; 



7T 

2 : 

K u sin 4>q cos i 



0. 



(7) 

(8) 
(9) 



The solution is found to be a Bloch domain wall of 
width A = \J Aj K u and energy a = AK u Nd w , where 
Ndw denotes the number of magnetic sites inside the 
domain wall. More precisely the domain wall profile 
is given by sinc/)o = sech [(x — xo)/A] and cos</>o = 
— tanh [(x — xo)/A]. 

To account for thermal fluctuations, we consider small 
deviations (69, 6<p) about the static configuration (9o, <j>o)- 
We expand H up to the second order with respect to 69 
and 6d> to obtain, 



6H = 



K u 



d 3 r {69 (V + k) 69 + 6<j>V5<j>} . (10) 



In agreement with earlier works J2SI2ZI we g nc j that the 
energy of the thermal fluctuations is described by a 
Schrodinger-like operator T> = —\ 2 d 2 — 2sech[(x — 
xo)/A] + 1 with k = K±/K u . The eigenvalues of V are 
and ojk = 1 + fc 2 A 2 . The zero-eigenvalue solution £i oc , 



Cloc(r) 



Ah-r 



V2N< 



sech 



dw 



x 



A 



(11) 



corresponds to the Goldstone mode of the system since 
the energy of the static wall is independent of its posi- 
tion Xo- In other words the £i oc -part of 8<f> contains no 
energy £i oc 2?£ioc = 0. We can avoid expanding 54> on 
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(x-x Q )/>. 



FIG. 2: Wavefunctions \/2Nd^^i oc (r) and y/i^>kN^k(r) about 
the domain wall et xq, for k x — 1/2A and k x — 4/A. 



the Goldtsone mode by elevating domain wall position 
to a dynamical collective coordinate XoM- The system 
is not rotationally invariant about the wire axis, because 
of the strong perpendicular anisotropy K± . As such, the 
£i oc -part of 89 carries a finite energy and the wavefunc- 
tion £i oc corresponds to a bound state of the system. In 
the following the amplitude of this bound state will be 
noted c\ oc (t). The non-zero eigenvalues of operator T> 
correspond to the propagating waves S,k(f), 



Mr) = 



1 



tanh 



X — Xq 

A 



(12) 



we have noted N the total number of magnetic sites in 
the sample. The wavefunctions are orthonormal to 
each other, 



d_ 3 r 
a 



3 CfcCm 



(13) 



which, in turn, are orthogonal to the bound state wave- 
function £i oc , 



d 3 r 

5~£fc6oc = o. 
a 3 



(14) 



The wavefunctions £i oc and are represented in Fig. [2] 

It is convenient to expand the small angle deviations 
89 and 8<p in terms of the eigenfunctions via the 
complex- valued variables dk through the transformation 
Sip + i59 = ici oc ^i c + J2k dk£k- Using this notation, we 
note that the Hamiltonian can be written as TL = a + 
W4) E fc + «/2) d* k d k - (k/4) (d k d. k + d* k d*_ k )}. 
This is not a quadratic Hamiltonian because of the finite 
perpendicular anisotropy Kj_, which leads to elliptical 
spin precession. To diagonalize this Hamiltonian, we 
follow the usual prescription by means of the Bogoliubov 
transformation c k 



l k d -k' with 



/ ujf, + k/2 ± Ml k /K u 

2m k /K u 



(15) 



where the frequency flk is defined as, 

2 



Wfe (uj k + k) . 



(16) 



Replacing dk and d* k by the magnon operators Ck and c k , 
we get the quadratic spin- wave Hamiltonian, 



SH = K ± c 



2 

loc 



k 



k c k c k 



(17) 



Spin wave energy is Mlk- The mode Ck has two compo- 
nents Ck — c k ci + c k h , which describe the wall deformation 
and the thermal propagating spin wave excitations. 

Next, we quantize the system by turning the com- 
plex variables c k h and c k h * into the boson operators Ck 
and c k , which obey the usual bosonic commutation rela- 
tions. For the sake of clarity, we will find it convenient 
to use the variables 4> k = c k + c*_ k , 9 k — (1/i) (c k — c*_ k ) 

and their corresponding operators 4>k — c k + c^ fc , 9 k — 
(X/i) ^Cfc — c^-kj ■ The small angle deviations 5<f> and 59 
are then expressed in terms of <pk and 9k as, 

84>{r) = (18) 

k 

S9(r) = c lo Mo C (r) + J2 e k4Ck(r). (19) 



The parameters v k = (u k +u k )/2 and v k = {u k —u k )/2 
represent the ellipticity of the spin precession. If the 
system were rotationally invariant about the wire direc- 
tion Kj_ = 0, then spin precession would be circular 

4 = * = 1/2- 

Lastly, it is convenient to renormalize the bound state 
amplitude c\ oc and introduce a new variable p as, 



V : 



A 



(20) 



We will see in the next section that p(t) represents the 
domain wall kinetic momentum. 



III. GENERALIZED ID MODEL OF BLOCH 
WALL DYNAMICS 

A. ODM without spin-waves 

As we will show in subsequent sections, the deforma- 
tion of the domain wall due to spin transfer torques is 
described by both the bound state c\ oc and the prop- 
agating states Ck- However in this subsection, we will 
disregard the spin wave modes Ck and discuss the role on 
the domain wall dynamics of the sole bound state am- 
plitude Qoc = p{t). We show that we recover the usual 
one-dimensional model of Bloch wall dynamics without 
spin- waves. 
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In the absence of the nonadiabatic spin-transfer term, 
it is possible to derive the equations of motion using a 
Lagrangian formalism. This is particularly well-adapted 
to the present problem in which the magnetic system 
is subject to constraints related to the presence of the 
domain wall. 

The total Lagrangian for the magnetic system is the 
difference between a "kinetic" (or Berry phase) term, 



Adn(u = 0) = S 



d 3 ? 



(1 — cos 9) d r 



(21) 



and the magnetic energy of the system Ti.. The inclusion 
of the adiabatic spin-transfer term appears as a moving 
reference frame at a velocity equal to the effective drift 
velocity of the spin current u. This is accounted for by 
replacing the time derivative in the kineti c part of the 
Lagrangian by a convective derivative] 39 * 42 ! 



-*kin 



s 



(1 — cos 9) (ud x + dt) 



(22) 



To zeroth order in the deformation, the kinetic part of 
the Lagrangian does not contribute to dynamics and can 
be neglected. The first-order term of £kin with respect 
to the deformation is 



(i) 



rf 3 7 



cioc6oc(m9 x + d t )(j>a, 



-p (u - ± ) 



(23) 



The overall magnetic energy is the sum of the static do- 
main wall energy a and of the dipolar energy K±c 2 oc . The 
latter can be re- interpreted as the kinetic energy p 2 /2m 
of the domain wall, 



By combining the Lagrangian functional £ (24 1 with the 
dissipation function T ( 26 ) , the Euler-Lagrange equa- 
tions lead to 



dp 
di 

dxo 

~dt ~ 



a2K± dxo 

aS dp 
'~ m2Kj_ ~dt 



= 0, 



P_ 
m 



(28) 
(29) 



Thus, the well-known equations of motion of the one- 
dimensional Bloch domain wall, Q and are recov- 
ered. We conclude that the bound state amplitude p(t) 
can be interpreted as the domain wall momentum as long 
as the propagating spin waves are neglected. 



B. Interaction between domain wall and 
propagating modes 

The expansion of the full Lagrangian C up to the first 
order in the spin waves ( 23 1 involves only the bound state 



and does not depend on the propagating waves. The in- 
teraction between the propagating waves and the domain 
wall arises from the second order expansion, 



/•(2) _ 
''-kin — 



d 3 , 



SS9 (ud x + d t ) 5<f>, 



\ ' 



1 



s jVkp-k 



^ Vkm9k<Pm 
k^m 



S(u - ± ) 



^ 2 

k 



C k Ck 



(30) 



,2 
loc 



P 

2m ' 



where m — Nd w S 2 / K±X 2 is the Doring mass. Finally 
the full lagrangian is obtained as 



C 



-p(u — xq) — a — p 2 /2m . 



(24) 



Gilbert or viscous damping can be accounted for by in- 
cluding the dissipation function, 



T = aS 



d 3 7 



) 2 + 6 2 sin 2 i 



(25) 



Assuming the "rigid" domain wall profile 9 = 9o(t) and 
<f> = 4>o(t), the dissipation function T is readily rewritten 
as 



:F = aS { 2 lPN dv 



(26) 



The equations of motion of the domain wall coordinates 
q = (xq,p) arc obtained by means of the Euler-Lagrange 
equations, 



dC _ d^dC _ \_dT_ 
dq dt dq 2 dq ' 



(27) 



The coefficient Vk describes a coupling between the 
bound state £o and the propagating waves, 



7T 1 & . ( kr.XTT 

Vk = ; vZ sech 

2 VUkN k 



(31) 



The coefficient Vkm describes a coupling between the dif- 
ferent propagating waves, 



_ e ^ itt (k x X) 2 - (m x X) 2 

Vkm — V k V m — - X 

Z L^JUk^m 

csch[ -Hk x + m x ) ' 



(32) 



with L = N x a the length of the wire. Coefficients Vk and 
Vkm are shown on Fig. [3] v k becomes negligible for k x > 
1/A, which means that the term — (it — io)^2 k Vkp4>k 
will mainly couple the domain wall to spin waves with 
a wavevector k x of the order of 1/A. In addition, the 
coefficient Vkm is large for long-wavelength spin waves 
fc < 1/A, to < 1/A and is roughly proportional to k x — m x . 
Thus the coupling term S(u - x )J2km v km9k(t>m will 
lead to a significant interaction with spin waves having 
k x ~ —m x ~ ±1/A. In some sense, the latter coupling 
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(a) 





term represents reflection of the propagating spin waves 
from the domain wall. 

The quadratic term S(u— ± n ) ^2 k (k x /2)c k c k represents 
a shift in the dispersion relation of the magnons. As such, 
the frequencies of the magnons depend on the relative 
velocity between the spin current and the domain wall, 



FIG. 3: (a) VNvk as a function of k x and k y . (b) Lvhm as a 
function of k x and m x ; K±/K u = 57. 




-1.5 -1 -0.5 0.5 1 1.5 

k k 



FIG. 4: Dispersion relation of the magnons as a function of 
the spin current velocity u, for K±/K u = 57 and h/K u = 1.88 
ns. 



£fe -> £fc - ttM" - &o). (33) 



This shift may be interpreted as a Doppler effect, which 
was originally put forward by Lederer and Mills^ and 
found recently in experiment by Vlaminck and Bailleul 44 . 
Some consequences of this Doppler shift have already 
been investigated so far, e.g. the excitation of mono- 
domains structures by a sole dc electrical currentP*! or as 
the enhancement of dissipation by spin transfer torque.^ 
An example of the current-induced Doppler effect on the 
spin-wave dispersion relations is shown in Fig. [4] We 
note that the current-driven terms can lead to a soft- 
mode with a nonzero wavevector k c , as indicated by the 
dashed line in Fig. [4] for u = 5A/ns for k c \ ~ 1.5. As 
pointed out by Shibata et al.p^this would lead to an in- 
stability in a uniformly magnetized ground state whereby 
the nucleation of domains of size ~ 7r/fc c is favored. 

Combining the kinetic part of the Lagrangian and the 
magnetic energy, the full Lagrangian C describing the 
magnetic system is found to be 



£ = Sy~] jO k (f>- k + S(u - x Q ) 



p 

S 



k 



s 



E 



h,-/- — rn 



t'kr 



2m 



J2 m kC* k c k . (34) 



The classical coordinates of the domain wall are coupled 
to the propagating spin waves through the interaction 
potential, 



V = S(u ~ xq) 



E 

L k 



Vk 



P 

s ( 



Vkr 



km 



(35) 



This potential must exist because the wavefunctions £o 
and £fc don't diagonalize the kinetic part of the lagrangian 
£kin- As such, a finite difference u ^ xq between the spin 
current velocity u and the domain wall velocity Xq will 
always give rise to a coupling between the propagating 
spin waves and the domain wall. However, if the spin- 
transfer from the adiabatic torque is achieved with an 
effective drift velocity such that u = Xq, then the interac- 
tion potential V will vanish identically and, in that case, 



the domain wall and the spin waves will be completely 
decoupled. We point out that the solution xq = u is ac- 
tually satisfied if the (3 coefficient in the ODM is identical 
to the Gilbert coefficient a, or in other words, the conser- 
vative and non-conservative dynamics of the system are 
invariant under a Galilean transformation. 



C. Force and torque 

In the following, we seek to generalize the ODM by the 
inclusion of the interacting potential ( 35 1 which couples 

The Eulcr- 



the spin-wave modes and the domain wall. 
Lagrange equations ( 28 1 and ( 29 1 become modified by 



the latter and therefore involve new terms with respect 
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to the magnons. 

By definition, the force F exerted on the domain wall 
is equal to the time derivative of the domain wall mo- 
mentum dp/ dt. The Euler-Lagrange equation 



dp 
~dt 



an _ d dv 

8xq dt dx a 



(36) 



indicates that F originates in both the magnetic energy 
Tl and the interaction potential V. 

Let us first consider the contribution Fa from the mag- 
netic energy Ji, 



F» 



on 



(37) 



This force may be rewritten as 
d 3 r 



-4K„ 



Sm z (x)d x sedx {x — x ) , (38) 



where Sm z represents the decrease in the longitudinal 
component of the magnetization Sm z — [SO 2 + S(/) 2 ) /2 
and is finite at finite temperatures due to thermal exi- 
ctations. By inspection of (38 1, we observe that Fa is 



only finite if 8m z is odd with respect to the domain wall 
position 5m z (xo — x±) = —Sm z (xQ + xi). This can be the 
case if the an electrical current flows through the domain 
wall and breaks the symmetry of the system. 

Let us now consider the force — 4i coming from the 



interaction potential V. 



dt dxQ 

This force may be divided into 



three different components: Fj, Fdef and -F s toc- The force 
F,, 



km 



(39) 



originates from the coupling between the different spin- 
wave modes induced by the electrical current. Because 
( 8k4>m ) depends on the statistics of the magnons, the 
force -F] is expected to depend on the temperature. The 
force -Fdef is related to the deformation modes c^ cf (see 
H. By writing ^ cf = c£ ef + c^f* , we find 



F 



dcf 



k 



(40) 



In contrast to Fj and Fd e f, the force F stoc is not de- 
terministic but stochastic. Its average value vanishes 
( F s toc ) = 0, but its autocorrelation function is finite 
( F stoc (t)F stoc (t') ) ^ 0. The force F stoc (t) may be rewrit- 
ten as, 



F s < 



d 
"dt 



v k p4> k 



(41) 



Adding all these forces together, the total force F exerted 
on the domain wall is found as 



dp 
dt 



Fu 



Fj + Fdef 



(42) 



It will be shown in the following that all the forces Fa, 
Fj and Fdef can be re-interpreted as a modification of the 
domain wall mass. 

So far we have looked at the contribution of the spin 
waves to the force F but still not to the torque T. The 
torque T is derived by taking the Euler-Lagrange equa- 
tions with respect to the momentum p, 



x = h u - 

m 



^2v k (f) k (x - u). 



(43) 



In addition to the spin-transfer drift velocity u and to 
p/m, which represents the torque due to the demagne- 
tizing field, the equation (43) involves a stochastic torque 
T 

± stoc 1 



(±0 - U) ^2 V k4>k 
k 



(44) 



This torque induced by the spin waves, which increases 
as a function of the relative velocity x — u. 



IV. SPIN- WAVE EMISSION 



A. Dissipation function 



According to (41 1 and (44 1, spin waves exert a non de- 



terministic force and torque on the domain walls. This is 
because the domain walls are sensitive to the fluctuations 
of the spin waves through the interaction potential 



(u - x )p ^2 v k (c k + cL 

k 



(45) 



We know by the fluctuation-dissipation theorem that 
these fluctuations necessarily lead to dissipation. As we 
haven't taken this dissipation into account yet, our ODM 
with the forces Fa, Fj, -Fdef, F'stoc and the torque T st0 c is 
not fully consistent. 

The spin current supplies energy to the domain wall 
motion by increasing the kinetic energy p 2 /2m. However 
a part of this energy is lost by the domain wall and is 
transfered to the lattice through a series of relaxation 
processes. The transfer of energy to the lattice accompa- 
gnying a domain wall motion is due to various channels 
of relaxation^! One channel of relaxation starts with the 
emission of magnons by the excited domain wall. 

Caldeira and Leggett,S3 following on from the semi- 
nal work and Feynman and Vernonj^ have shown irre- 
versibility to arise when a moving particle is coupled to 
numerous degrees of freedom. In our case, the magnetic 
domain wall represents the moving particle, which prop- 
agates through a dissipative environment represented by 
the thermal spin-waves. The domain wall is described by 
its position xq and its kinetic momentum p, whereas the 
spin wave environment is described by a set of oscillators, 



Wbath — ^ fiil k c k c k . 

k 



(46) 
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The bath of magnons is assumed to be always at ther- 
mal equilibrium which is reached through very fast three 
and four magnon processes and due to the interaction 
with phonons. The domain wall and the spin wave envi- 
ronment are coupled to each other through the "linear" 
interaction potential Vi n t. Notice the similarity between 
the interaction Vj nt and the interaction between electrons 
and phonons in the metals. The present system is ac- 
tually formally equivalent to the dissipative system in 
Caldeira and Leggett's general theory which itself is in 
agreement with the fluctuation dissipation theorem! 49 ! 50 ! 
Notice that a coupling very close to Vint has already been 
considered by Thomson and Stamp, who investigated the 
damped motion of vortices driven by a magnetic field 
The emission of a magnon k is represented by the term, 



V e m = (u-X )pV k cl. 



(47) 



Similarly the absorption by the domain wall of a magnon 
k is represented by, 



(48) 



The energy transfer from the domain wall to the spin 
waves ensemble will correspond to the difference between 
the emission of magnons and the absorption of magnons. 

In the following, the various states of the spin wave 
ensemble will be denoted by \E n ) and the density of 
state of the spin waves will be represented by p(E). By 
Fermi's golden rule, if the spin waves ensemble were in 
the state \E n ), then the exchange of a magnon k be- 
tween the domain wall and the bath of the spin waves 
would statistically decrease the domain wall's energy by 
an amount F n (k) as, 



F n (k) = 2nn k 



Ei 



Ern\V em \E„ 



2 s(E m - e„ - nn k ) 



- I ( E m \V* b \E n ) \ 2 8{E m - E n + m h ) 

m 

which can be written explicitly as, 



F n (k) = 2nSl k 



\(u - x )pv k \ 2 S(E m - E„ - e k ) 



J~] \{u - x )pv k \ 2 5(E m - E n + e k ) 



.(49) 



In reality the probability for the spin wave ensemble to 
remain in the state ( E n ) is equal to the Boltzmann fac- 
tor P n = exp (—(3E n ). The magnons k thus contribute 
to a transfer of energy from the domain wall to the spin 
waves ensemble by an amount F(k) — ^2 n PnF n (k). In 
the continuum limit Y2 n E n — > J dE, 



F(k) 



dEp(E)f(E)F E (k) 



(50) 



where f(E) = exp(—/3E) and F E (k) is, 

F E (k) = 2-Kfl k \(u ~ x )pv k \ 2 x 

[S(Em -E-m h )- 8(E m -E + hfl k ) 



In other words the energy lost by the domain wall 
through the transfer of a magnon k is, 



F(k) = \(u-xo)p\ 2 nlv 2 k R(k) 



(51) 



where the function R(k) is closely related to the spectral 
function of the environment and equals 

POO n 

R{k)= / dEp(E)f(E) — [p{E + e k )-p(E-e h )}, 
Jo "fe 

with p(E) being the density of states p(E) = J2 m S(E m — 
E) of the spin wave ensemble. The total energy lost by 
the domain wall per unit of time J- = J^ fc F k is finally 
obtained as, 



F=J2\(u~io)p\ 2 n 2 k v 2 k R(k). 



(52) 



Let us introduce the dimensionless friction paramter 
i] = J2 k ^k v kR(k)- The dissipation function can then 
be rewritten as, 



T 



2*7 
N dw S 



\(u-x )p\ 2 . (53) 



B. f3 term 

To understand how this dissipation function affects 
current-induced domain wall motion, it is convenient 
to perform a Galilean transformation and use the local 
frame moving at the velocity of the spin current. The 
domain wall position in this frame is 



q = xq — ut. 



(54) 



The dissipation function T = —dTL/dt can be rewritten 
in terms of q as, 

thus treating q and p as independent variables, 

T= 2-dq q (56) 

If the system is translationally invariant (no extrinsic pin- 
ning), the lagrangian describing the domain wall motion 
will not depend on the the domain wall position q 



£ = -H+—q 
oq 



(57) 
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therefore the time derivative of the domain wall Hamil- 
tonian becomes, 



dH 

lit 



q- 



(58) 



d dC _ dC\ 
K dt dq dq J 

Combining ( 58 ) with the definition of the dissipation 
function T ', we get 

d^dC 
dt dq 



dC 
dq 



1 dT 



2 dq 



= 0, 



(59) 



which after carrying out (24) and (55 1 yields, 

. ^ 2rip 2 . 2rip 2 
NdwS Nd w S 



(60) 



The dissipation by spin wave emission gives rise to a 
damping term 2np 2 xo/Nd w proportional to the domain 
wall velocity ±0 and to a force 2r)p u/Nd w S proportional 
to the spin current u. These are equivalent to a damping 
coefficient a sw and to a coefficient (3 SW expressed by, 



0s 



Tjp" 



N dw K±m 



(61) 



It is seen that the spin wave contributions a sw and (3 SW 
are identical to each other and are also proportional to 
the domain wall kinetic energy p 2 /2m. 

The dissipation function generally tends to decrease 
during the motion, because it vanishes at equilibrium. In 
the case of Gilbert damping, which is described by the 
dissipation function F/2SNd w = a(io / X) 2 (1 + a 2 ) , dissi- 
pation thus tends to reduce the domain wall velocity xq 
to zero. In contrast, the spin waves which lead to the dis- 
sipation function ( 55 1 will not tend to lower the domain 
velocity Xq, but instead the relative velocity q = \x — u\. 
In other words, whereas Gilbert damping slows down the 
mowing domain walls and reduce their velocity to zero, 
the damping by the spin waves acts towards restoring the 
solution Xq = u. 

The equality between a and (3 restores Galilean invari- 
ance, which has previously been argued by Barnes and 
MaekawaP^I In our theory, this invariance is restored be- 
cause the dissipation channel, in this case the magnons, 
also "flows" with the effective drift velocity u through 
the Doppler effect. As such, the dissipation of the do- 
main wall motion through spin-wave interactions leads 
to a Landau-type damping in which the motion tends 
towards xq — u. We contend, therefore, that any dis- 
sipation channel that drifts at the same velocity as the 
spin-current u should lead to a similar dissipation chan- 
nel, for which we can write symbolically as adrift = Airift- 

Of couse, extrinsic pinning centers such as magnetic 
impurities break the translational invariance of the wire, 
so ( 60 1 would no longer apply and the full coefficients 
a and (3 should not be equal in general. Notice that 
the damping coefficient a sw , which we have identified in 



(61 1, doesn't exactly recover the Gilbert coefficient a, 
since it doesn't appear in the Euler Lagrange equation 
with respect to p. 
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FIG. 5: (a) Correlation function I(t—t') between the stochas- 
tic torque T st0 c(t) at two different times t and t' , for a wire 
of Permallo^Pl, with K±/K u = 57 and h/K u = 1.88 ns. (b) 
I(t - t') for a wire of C0/P1P, with K ± /K u = 0.73 and 
h/K u = 0.013 ns. 



Stochastic field 



The force F st0 c (411 and the torque T st oc ( |44| , which 
accompany the dissipation by the magnons, lead to a do- 
main wall motion that is stochastic. As ( <f>k ) = and 
( 0k ) = 0, the average values of F stoc (t) and T stoc (t) van- 
ish identically, but their correlations are finite. Assuming 
the wire to be very large compared to the domain wall 
width A, the auto-correlation of the force F stoc is found 
to be 

( w*)*itoc(0 > = pWp(0 (^) 3 ^/(* - 



■ P (tW) 



Ku 

h 



( j L _ V 2k B T m _, 



V27TA/ Ku 



with, 



I{t-t') = 



d 3 k\ 3 (4f^ 



j(t-f), 



k T \ir 



sech 



x CO f^-^ , (62) 



Wfe (uJk + k) : 



and 



j( t -t')= [d s k\ 3 (vif(iy 



sech 



k x \TT 



cosfi fe (*-i')- (63) 



I(t — t') is presented on Fig. [5 with the anisotropy con- 
stants K u and K± taken to be those for permallo}^ 
and Co/PtPM It is seen that I(t — t') cancels out very 
quickly when \t — t'\ becomes larger than 1 ns. Be- 
sides the function J(t — t') behaves very similarly and 
weakens very quickly when the time difference \t — t'\ in- 
creases. We therefore conclude that the stochastic forces 
-Pstoc(i) and F stoc (t') are not correlated to each other at 
the characteristic time scale associated with domain wall 
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motion, which usually is longer than the nanosecond.^ 
For instance, the timescale of a domain wall at a typi- 
cal velocity of 1 m/s in a permalloy wire corresponds to 
X/lm.s^ 1 ~ 50 ns. 

Analogously, the auto-correlation of the torque T stoc (t) 
is found to be 



( -^stoc(^)^stoi 



l(0) = MfJL) J » (t . 0) 



V2ttA/ K u 



and the inter-correlation between the force F stoc (t) and 
the torque T s t ac (t') is found to be 



to \27rA/ A M 

These correlations involve the same function I(t — t') as 
before (62 1. From this study on the fluctuations of the 
spin waves, we conclude that the stochastic force -Fstoc 
and the stochastic torque Tgtoc behave as a multiplicative 
white noise on the domain wall. 



> 3 2k B T . 



V. REDUCTION OF DOMAIN WALL WIDTH 

External forces, such as a magnetic field or an elec- 
trical current, lead to deformations of the domain wall 
profile. The most evident feature of such deformations is 
the creation of dipolar charges, which collectively build 
up the kinetic energy of the domain wall. Indeed the ki- 
netic energy K±X 2 p 2 /2Nd w S 2 depends on the transverse 
anisotropy K±, which in permalloy is due to the surface 
dipolar charges on the thin film surface. Another feature 
is that the width of the domain wall is reduced by the 
external magnetic field throughout its propagation^. In 
this section it is demonstrated that the latter feature does 
not specifically apply to the drag by a magnetic field but 
also appears in a similar manner under applied electrical 
currents. 

The domain wall deformation is entirely de scribe d by 
the modes cj oc (t), c k (t) and <4(t). In Section ( III A I the 
mode c\ oc (t) has been interpeted as the kinetic energy of 
the domain wall. In the present section, we show that 
the modes c^ of (i) and cjf c! (t) can be interpreted as a 
reduction in the domain wall width, which subsequently 
appears as a change in the domain wall mass. 

Departing from the Lagrangian ( 34 1 and considering 



the Euler-Lagrange equation with respect to ct, we ob- 
tain 



5 
2i 



Cfe 



- (hQ k - S(u- ± )y^) Cfe - (u - x )v k P 
-{u-± Q )S ^2 [ 



k^m 



^m, — k (c?n C_ rn j 



(64) 



The mode c k has two components c^J 1 and c^ cf . The fast 
component ^(t) (~ 4 GHz in permalloy) represents the 



thermal excitations, whereas the slow component c^ cf (~ 
20 MHz for dx^j 'dt=lvs\/s in permalloy) represents a de- 
formation. As the differential equation ( |6~4| is linear, the 
components c fe h and c^ of can be calculated independently 
The slow component c£ is almost static c£ cf = and 



corresponds to the "static" particular solution of (64 1 



In contrast the fast component of c k has no static part 



Let us calculate the deformation mode c2 from (64 1 



and is an homogeneous solution of ( 64 1 
he < 

by approximating c k de ' 
Taylor 6 
locity u — in 



= 0. It is useful to consider the 
Taylor expansion of c^ ef with respect to the relative ve- 



dcf 



! n (u - XqY 



(65) 



Since the domain wall is much slower than the propagat- 
ing spin waves (u — Xq)/X <C Q k , it is sufficient to keep 
only the first order in the Taylor expansion ( 65 1 , 



c fc ef »0i(u-£o). 



On readily obtains 



„def 



(u - x ) v k 
Ml k 



-P- 



(66) 



(67) 



Therefore the deviation S(f) dcl in the spherical angle <fi due 
to deformation is 



/def 



E K { + c -t 



/def 



= -*E 



x sech 



(u - x a )p 
Ml k 



1 



(68) 



(69) 



By letting (uf\ = (l/4)y/(uj k + n)/u; k and p 
(±o — u) to, one finds 

7T ^ S 2 (u-x ) 2 1 



idef 



4 ^ 

k 



x ^sech 
a 



K U K ± X io k 
k r Xir x 



1 



a ■ (70) 



2 ) yfiJk~N 

Let us now calculate the deformation angle 5(j) de f 



cor- 



responding to a change SX in the domain wall width. Its 
comparison to ( |70| will then allow us to estimate SX. 
The (^-angle describing the magnetization inside a Bloch 
domain wall is. 



</>o(z) 



cos 



tanh 



x 



An 



(71) 



If the domain wall width is modified by the spin current 
by 5X(u), the actual magnetization profile <j) will deviate 
from the equilibrium Bloch profile </>n as 



x — Xo 
~~X 2 ~ 



lanh.2 f ] §X , (72) 
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The deformation angle 5<fi de f {r — rn) = <p(r) — 4>o(r) can 
be expanded on the wave functions £j.(r — t*o) as 



ro 



^2Pk£,k(r 



(73) 



where 
Pk = 



Scf>(r)C-k(r) 




tanh 



x - x 
An 



Only the real part of the wave function £_fc(r) will con- 
tribute to Pf. , 



Re[£_ h (r)] = 

accordingly, 
Pk 



1 



cos [k x {x — xq)] tanh 



x — Xq 

A 



+fca;Asin fe x (x — xq)] cos(fcj^.r), 



c5A iVi7r /" fc x A7r 

: sech 



(74) 



Hence the change 5 A in the domain wall width corre- 
sponds to the deformation angle, 



/dcf 



ESX Nj_Tr ( k x \Ti\ 



(75) 



Comparing ( 70 I and ( 75 1 , and approximating uj^ m 1 , we 
finally obtain 



SX 

T 



lS 2 (u-x ) 2 
"4 F M F ± A 2 



(76) 



Eq. ( 76 1 shows that the domain wall shrinks under the 



effect of spin transfer torque. According to (40), this 



reduction in width leads to a force Fj c f which is non- 
linear with respect to the relative velocity u — in, 



F 



clef 



-c— 
dt 



S 2 {u-x ) 



2 1 



P- 



(77) 



with c = (tt/32) /_ ~ dk x sech 2 (Trfc^/2) /(l + k 2 x ) ~ 0.1. 
In the case of permalloy nanowires, represents the 
demagnetizing energy and cancels with the factor S 2 in 
the numerator in ( 77 1 . Therefore F^ e f does not depend 



much on the transverse anisotropy K±. 

The deformation force -Fdef can be re-interpreted as a 
modification of the kinetic momentum and of the mass of 
the domain wall. As the domain wall deforms, its kinetic 
momentum becomes 



p + Sp = p[ 1 



3c5* 2 (w-x ) 2 
K±A 



and its mass is modified by 
8m 6c 



A^ 



(u - x y 



9VB 



(78) 



(79) 
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FIG. 6: Modification of the domain wall mass and of the 
domain wall width as a function of the relative velocity \u— xo\ 
and the exchange stiffness constant A. 



The modifications of the domain wall mass and of the 
domain wall width are represented on Fig (6] as a function 
of the relative velocity \u — xq\ and for different values 
of the exchange stiffness constant A. They are almost 
identical 6m/m ~ —5X/X. 



VI. TEMPERATURE DEPENDENCE OF THE 
DOMAIN WALL MASS 



The flow of an electrical current through a ferromag- 
netic wire modifies the spectrum of the magnons accord- 
ing to ( 33 ) and thereby shifts their average kinetic mo- 



mentum to a finite value k = ^2 k kriB(k)/ ^2 k Jig(k). 
The applied electrical current therefore leads to magnon 
current. Inside a domain wall, magnons behave slightly 
differently because of the sing ular coupling potential 
—S(u - Xo) J2k m v kmQk4>m <|34|). This potential may be 



interpreted as a modification of the domain wall's en- 
ergy due to its interaction with the spin wave environ- 
ment. This section is devoted to this coupling potential, 
whereby we investigate its consequences on the domain 
wall dynamics. Specifically, this potential will be shown 
to contribute to the forces Fj ( [39] ) and Fh (38) on the 
domain wall. Some light will thenbe shed on the process 
suggested in the Introduction (Section [I]) . 

The relationship between the forces Fj and Fh and 
the statistics of the magnons will be established in Sec- 
tion VI A The force Fj arising from the dc-component of 
u will be carried out in Section [VI B| and the force Fh aris- 
ing from the ac-component u(t) will then be calculated in 
Section VI C Both these forces will be re-interpreted as 
a modification of the domain wall's effective mass. The 
domain wall mass will in turn become sensitive to the 
actual temperature of the ferromagnetic wire. 
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A. Forces Fj and Fh 

While the force Fj depends on the correlation 
between 89 and 8(f), the force Fh depends on the auto- 
correlation ( 89 2 ) and ( 5(j) 2 ). These correlations, which 
are finite under the effect of an electrical current, can be 
calculated with perturbation or linear response theory. 
For the sake of simplifying the notation, we will employ 
in this section a 2 x 2 matrix representation which is de- 
tailed in the Appendix. 

The correlations between the angular deviations 8<f> 
and 89 are well described by the lesser magnon Green 
function T) < (r,t,r',t') defined as 



ihV < 



69(r',t')69(r,t) 
56(r',t')6(f>(r,t) 



(r',t')66{r,t) 
(r',t')S<f>(r,t) 



The representation of the lesser function in the momen- 
tum space T> < (k, t, k' , t') is obtained by expanding 69 and 



8(f) on the wavefunctions according to (18) and (19 1. 
The representation in momentum space is quite useful as 
it reveals the relationship between the angular deviations 
8(j>,89 and the statistics of the magnons. 

The force Fu depends on the diagonal component of 
the lesser function = Vf a = T>^,, 



F K = J2f H (k,q)ihV<(k^,t,k +1 t), (80) 

hq 

where k- = k — q/2 and k 4- q/2 and 



1 + 3(k,\f + 1^ 



u u i U> U> 



The force Fj depends on the off-diagonal component of 
the lesser function 2?< ff = 



Fj=J2f( k ^)j t ^ s (k-,t,k + ,t), 

hq 



with 



-Sv-fc_,fc 4 



(82) 



(83) 



We recall that the coefficient was defined in (32) 

and that the coefficients Vki v< t were defined in (IIIBl. 

As / H (fc,q) is odd in q and even in k, the force Fh 
will be finite if Djp (fe_ , t, k + ,t) is odd in q and even in k. 
It will be shown in Section IVI CI that the non-adiabatic 
response of the spin waves to a dynamical spin current 



u(t) gives rise to such a lesser Green function T> < . Besides 
as f ' 3 (k, q) is odd with k and even with q, Fj will be finite 
if T>'^ ff (k_ 1 t 1 k + ,t) is odd with k and even with q. We 
will see in ( VIB[ ) that this is the case for the adiabatic 
response of the spin waves to a dc spin current u. 
k- |,w 



Awvv 



2' L 



FIG. 7: Representation of the force Fj created by the spin 
waves when the relative velocity u — xq is finite and the mo- 
tion is adiabatic. The wavy line represents the magnon q 
exchanged with the domain wall, the vertex x arises from the 
force Fj and the vertex • represents the potential V. 



B. Adiabatic response 

The process in Fig. [7] depicts the adiabatic transfer 
of kinetic momentum between the propagating magnons 
and the domain wall, which arises from the coupling po- 
tential V ( 35 ) . Due to adiabaticity, the energies of the 



incoming magnon fe — f and the outgoing magnon k + f 
are identical. 

The statistics of the magnons perturbed by the cou- 
pling potential V can be investigated by expanding the 
contour ordered magnon propagator T>(k— f , k+ f , r, r') 
up to the first order in V, 



2?W(fc_,T,fc+,T / ) = -S(u-i ) [ dr x 

Jc 



w-fe_,fe + (n) 



^°\k_,T-T 1 )cj y V^{k + ,T X -T , ).{M) 



In ( 84 1 the superscript (0) represents equilibrium (i = 
u), whereas the superscript (1) represents first or- 
der perturbation theory. The lesser Green function 
2?(fc_, k+,t — t') can be obtained from the time ordered 
Green function with the Langreth formula (I?I?) < — 
T> r T> < + T> <r D a and the fluctuation dissipation theorem, 



, lo) = n B (uj) \ T>^ r (k, lo) ~ V^ a {k, oj)J . (85) 
By noting Im(z) the imaginary part of z, and 











i 









I = 



the lesser Green function is found in terms of the retarded 



(A7l and advanced (A8l propagators as 
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V^<(k-,k+,u) = 45(u-io)«-fc_,fc + n B (w) 



/Im^f r (fe_, W ) 5 ^(fe + ,c,)- g ^(fc_, W ) 5 ^(fe + ,c,) 
+ a„Im (gf r (fe_ , u)gf r {k + , u) + g^ r (k_ , u)g^\k + , u)) 



joy 



(0)r, 



(86) 



From this equation, we can infer that 
I?J| 1 ' )< (fc_, t, k +1 1) = 0. Thus by inspection the 

adiabatic response of the spin waves to (it — xq) doesn't 
contribute to the force Fh ( 80 1 . 

diagonal component T> i ^ < {k_,k^ 



In contrast the off- 
, t) of the lesser Green 



function is finite when the domain wall velocity differs 
from the velocity of the spin current xq ^ u. It is 
obtained as 



ifiV^ < (k-,k + ,t) = -8nS(u - x )v. 



k- ,fe_j. 



dni 



(<*) 



and leads to a force Fj ( 82 ) 
dp 8ttK^X 2 ^ 



F 



dt N, 



dw 



k.q 



dn B (e k ) , ^2 
— ^ — («-*_,*+) ■ 



(87) 



Being proportional to dp/dt, the force Fj may be re- 
interpreted as a modification of the domain wall's mass. 
In thin films, where the magnetization is uniform over 
the thickness of the wire, the change in the domain wall 
mass is found to be 



Sm ksTa 2 2ir 
m A N z 



where N z represents the number of atomic layers within 
the wire's thickness, a is the lattice constant, L y stands 
for the width of the wire and ^{K,L y ) is a function of 
L y and the anisotropics k = Kj_/K u . The change in the 
mass Sm /m is presented in Fig. [8] as a function of the 
wire thickness h — N z a and for different values of the 
temperature. In a permalloy wire with a cross section 
2 nm x 100 nm, the mass of the domain wall varies by 
- 5% between T = K and T = 600 K. This variation 
of the domain wall mass could be much larger near the 
Curie temperature T c ~ 750 K. 



C. Non-adiabatic response 

The process in Fig. [9] shows the non-adiabatic transfer 
of kinetic momentum between the propagating magnons 
and the domain wall, which arises from the coupling 



potential V (35) when the domain wall is accelerating. 



Through this process, the energy of magnons is not con- 
served. 

In the following, the domain wall motion is assumed to 
be harmonic u(t) oc e . This assumption will allow us 




h (nm) 



FIG. 8: Modification of the domain wall's mass Sm/m as func- 
tion of the wire's height h and the temperature T, assuming 
57. 



2> L 



n 

2 




2 ,U! + 2 



FIG. 9: Representation of the force Fh created by the spin 
waves when the excitation of the spin waves by the spin cur- 
rent is not adiabatic. The wavy line represents the momentum 
q exchanged with the domain wall and the energy Q trans- 
fered by the spin current, the vertex x arises from the force 
Fh and the vertex • represents the mixing potential V. 



to calculate the linear response of the spin waves to the 
domain wall acceleration in a more tractable way. The 
frequency of the domain wall acceleration (e.g. ~ 10 
MHz) is much lower than the spin waves eigenfrequencies 
fife (~ 5 GHz), but still induces an overlap between the 
different spin waves modes. As the actual overlap be- 
tween the spin waves depend on their spectral linewidth, 
the lifetime of the magnons is a key physcial quan- 
tity in this problem. According to the phenomenologi- 
cal Gilbert damping, each magnon k has a lifetime 
inversely proportional to its eigenfrequency and to the 
Gilbert coefficient, t^ 1 = aflk- The linewidth of each 
spin wave mode is therefore about A/ = 0.01 x 4 GHz ~ 
40 MHz. As the wire length L x is generally much larger 
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than the domain wall width A, the energy gap between 
the peaks of two consecutive spin wave modes is of the 
order of y /K^KZ(2irX/L lc ) 2 /2 ~ 5 MHz. 

For the sake of determining the linear response of the 



spin waves to the domain wall acceleration, we use the 
Dyson equation (84) as a starting point. As il/^fc <C 1, 
we linearize pW^with respect to fl, 



p(i)<(fe-|,fe + |, w ,o; + 0)=-m_ fe+fife+f (-n) 



dn B 



(89) 
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FIG. 10: Contribution of the transverse spin wave modes n y 
to the change in the domain wall's mass 5m /m. L y — 400 
nm, a = 0.01 and k = K±/K u = 57. 



The diagonal component of the lesser magnon Green 
function is then obtained as, 



D | ( | 1)< (fe_,i,fc+,t) 



with 



*(*..«-) = ZM.AcBch 



aS dp 
m dt 



n\q x 



k B T<&(k Xl q x ) 



(90) 



47T 



duj 



((hu - e k _) 2 + (ae fc _) 2 ) JJfiu) - e k+ ) 2 + (ae k+ ) 2 ) ' 
The force Fjj due to this non-adiabatic interaction is de- 



rived by carrying over the magnon Green function ( 90 1 to 



the definition of the force ( 80 1 . Because it is proportional 



to dp/dt, the force Fh can be re-interpreted as a modi- 
fication of the domain wall's mass. The contributions of 
the transverse spin wave modes n y to 5m/ m are shown 
on Fig. [lO] for different values of the temperature. The 
latter figure indicates that the change in the domain wall 
mass is in general less than 1/1000. This implies that the 
non-adiabatic response of the spin waves to the domain 
wall acceleration does not affect domain wall motion sig- 
nificantly. This is in contrast with the adiabatic response 



VII. DISCUSSION AND CONCLUSION 

In conclusion, we have developed a theory of current- 
driven domain wall dynamics which includes spin-wave 
interactions. In general, a spin current u traversing the 
wire will induce an interaction between the domain wall 
and the propagating spin wave modes. This interaction 
is proportional to the kinetic energy of the domain wall, 
m(u — ±o) 2 /2. As this kinetic energy is governed by the 
difference between the spin current u and the domain wall 
velocity xq , a domain wall moving at the spin current ve- 
locity will not interact with the spin waves. This coupling 
between the domain wall and the propagating magnons 
is shown to have two main effects on the domain wall 
dynamics: 1) Damping of the domain wall motion; 2) 
Renormalization of the domain wall effective mass. The 
damping accompanying the emission of magnons tends 
to decrease the energy of the domain wall thus tends to 
restore the solution, where the domain wall velocity and 
the spin current are identical xq = u. This dissipation 
channel is analogous to Landau dissipation in plasmas. 
The renormalization of the domain wall effective mass 
has two origins: the reduction in the domain wall width 
and the renormalization of the domain wall energy by 
thermal magnons. 

In the presence of Gilbert damping and additional 
spin- wave interactions, our extended one-dimensional 
model of Bloch wall dynamics is given by 



2K ± m* 
S 

2K±m 
S 

_P_ 

m ' 



VP 



N dw K±m* 



x 



^ + N dw K±m* U 
aS 

Xq — u - 



2m* 



(91) 
(92) 



of the spin waves ( VI B ) 



The effective mass m* takes into account the deformation 
of the domain wall and its coupling to the magnons. As 
such, it depends on the actual temperature of the wire. 
The emission of magnons, which is parametrized by the 
dimensionless coefficient 77, is important if the coefficient 
1] is close to 1 . The equation of motion ( 92 1 indicates that 
the emission of magnons "contributes" to the (3 term. 
Since in ( 92 1 the Gilbert damping coefficient a becomes 
shifted as a — > a+r]p /Nd w K±m* , emission of spin waves 
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emission of spin waves 
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FIG. 11: Proposed experiment for spin wave emission. A 
dc power supply Vd c generates a charge current j c , which in 
turn excites a domain wall pinned inside a wire. As the spin 
current u is finite but at the same time the domain wall is 
fixed, the relative velocity u — Xo is finite and spin waves are 
emitted. 



also "contributes" to the damping. In the limit in which 
the only dissipation channel for the domain-wall walls is 
through spin-wave interactions, we find that the Galilean 
invariance in restored and a — f3. 



The equations ( 92 1 and ( 92 1 describe current induced 



domain wall motion by including some mechanisms that 
depend on the temperature. Temperature takes effect 
on domain wall dynamics through the domain wall mass 
m* and also possibly through the emission coefficient rj. 
According to Yamaguchi et alP^, temperature inside the 
ferromagnetic nanowires may approach the Curie tem- 
perature T c because of Joule heating. According to our 
theory, such an increase in the temperature would appear 
as a renormalization of the domain wall mass. Our calcu- 
lation has indicated a change in the domain wall mass of 
about 5% between T = 600K and T = OK. This change 
becomes even more important close to the Curie temper- 
ature T c . Laufenberg et alP^l found a significant decrease 
in the spin transfer efficiency, when the temperature was 
increased by few hundreds of Kelvin. These authors have 
suggested emission of magnons to be responsible for this 
loss of efficiency. According to our theory, emission of 
magnons would instead help to depin the domain walls. 
Thus the experimental results of Laufenberg et al. can- 
not be explained within our model. 

According to the present theory, current-driven do- 
main walls might be viewed as the generators of spin 
waves. One experimental realization of such spin- wave 
generation is proposed in Fig. |11| where we show a do- 
main wall pinned inside a wire that is traversed by a dc 
spin current. The action of the spin-current, as we have 
discussed in some detail above, acts to deform the do- 
main wall and give rise to a non-equilibrium magnetic 
configuration. As its velocity vanishes in = 0, its energy 
m*(u — ±n) 2 /2 will become very large. A part of this en- 
ergy will then be dissipated to the medium through the 
emission of magnons. 



APPENDIX A: PROPAGATORS 

The magnon propagators allow the linear response 
of the spin waves to be derived in a very tractable 
way. The magnon propagators are defined by means 
of the magnon operators 4> k = c k + cl fe and k = 
j [ck ~ $-kj ■ Specifically the time ordered magnon 

propagator T> a p{k,k' ,t,t') for a, {3 £ {0,(j>}, is defined 
as 



T> aP (k,k',t,t') 



with 



Ta k (t)0-v(lf) = 



■j_(\Ta h (t)/3_ k ,(t') |) , (Al) 



0-k>(t')a k (t) if t < t' 



It is useful to let the times t and t' belong to a contour 
in the complex plane, which allows the Green's function 



( Al I to contain information on the non equilibrium prop- 



erties of the spin waves. The lesser and greater propaga- 
tors are defined as 

■D<p(k,k',t,t') = - p- k ,{t')a k {t) |), 

■D>p(k,k',t,t') = - l j_{\a k {t)p_ k ,{t') |), 

and the retarded/advanced propagators are 

v^(t,t') = e(t-t')(v> (t,t')-v< p (t,t>)), 
v a a0 (t,t>) = 9(t'-t)(v< (t,t')-v> (t,t')). 

To simplify the notation, we represent the various 
magnon propagators by means of 2 x 2 matrices as 

T>(h h 1 t t'\ - ( v s,e{k,k' ,t,t') V e ,4,(k,k',T,T') 
1 ' y V^ e (k, k', r, t>) V^(k, k', r, t') 

We use the interaction picture with respect to the unper- 
turbed spin-wave Hamiltonian TCq = X)fc^fc c l c fc- The 
creation and annihilation operators of the magnons have 
the time dependence c k (t) = c k e~ tnht , c k (t) = c* k e %nht . 

Before calculating the response of the spin waves to 
the coupling potential V, the unperturbed propagators 
need to be determined. The free retarded and the free 
advanced magnon propagators are obtained as 



V^ r (k,L0) 



I - er„ 



and 



V^ a (k,L0) = 



huj — e k + iO hui + 6fc + iO 



I + (T v 



huj — e k — iO huj + e k — iO 



(A2) 



(A3) 



The free lesser Green function is given by the fluctuation- 
dissipation theorem, 

vW < {k,uf) = n B {u) (V 0)r (fc,w) -£>(°) Q (fc,w)) .(A4) 



16 



(a) 



(b) 




FIG. 12: (a) Odd part (vkm — «mfc)/2 of coefficient vu m - (b) 
Even part (vkm + ^mfe)/2 of coefficient Vhm- 



It turns out to be useful to introduce the electron-like 
propagators g^ r (k,u)) and g±^ a (k,u}), denned as 



fiw =F e fc + ihTl 



1 



(A5) 
(A6) 



These propagators are related to the magnon propagators 
dA2| and (|A3| through 



gf a {kM = (gf r {k 



V** = (I + a y )gf r (I a y )g 



,(0)r 



and 



= (I+a y )gf a -(I-a y )g« 



(0)a 



(A7) 



(A8) 



The products of the retarded and advanced free magnon 
propagators, which appear throughout the calculation 
of the response, can be readily computed by using 
fA8|,([A7|, and by noticing that (7 + a y ){I - <r y ) = 



X(I±a y )i=2(I±a y ). 
The scattering potential V p5|) may be divided into 



a "linear" term and a "second order" term with re- 
spect to the propagating modes. The "linear" term 
does not perturb the correlations of the spin waves since 
( TcLck'Ck' ) = and ( Tc^Cfe/cL ) = 0. However the 
second order term —S(u — Xq) J2km v km^m4>k modifies 
these correlations. In the following we assume the coeffi- 
cient Vk m to be odd, 



Strictly speaking this assertion is wrong in general when 



k ^ 0, however the plots in Fig. 12 suggest that this is 
a good approximation. In the interaction picture, the 
contour-ordered magnon propagator is given by 

V a p{k,k',T-T') = -^(\Texp(-~^dT 1 V(T 1 )) 

x a fe (T)/3_ fe ,(r') |) . (A10) 

Linear response theory focuses on the first order expan- 
sion of T> a p(k, k' ,r — t') with respect to the coupling 
potential V. We use the superscript W to denote the 
first order expansion of the magnon propagator. Recall- 
ing Wick's theorem and using the relationship (A9|, we 
readily obtain 

pW(fc,fc',r,r')= / rfr lU _ fefe ,(r 1 )(pi° ) (fe,T-r 1 ) 

xv$(k', n - t') - 2>£>(fc, r - n)V$(k', n t>) 

Employing the 2x2 matrix notation, we finally find 

pW(fe,fe',r,r') = / dniu-kk'in) (All) 
Jc 

xVW(k,T-n)a y vW(k',n-T'). (A12) 



Equation ( |A12 1 expresses the linear response of the spin 
waves to the coupling potential V ( 35 1 . 



Vkr 



(A9) 
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